Local electromigration model for crystal surfaces 
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We analyze the dynamics of crystal surfaces in the presence of electromigration. From a phase 
field model with a migration force which depends on the local geometry, we derive a step model 
with additional contributions in the kinetic boundary conditions. These contributions trigger various 
surface instabilities, such as step meandering, bunching and pairing on vicinal surfaces. Experiments 
are discussed. 



In the presence of an electric current, mobile atoms ex- 
perience a diffusion bias, which is called electromigration. 
Surface electromigration is known as major source of rup- 
ture of microsctructures, and is also a interesting tool 
for spontaneous nano-structure formation. In the past 
15 years, kinetic instabilities at the micro, or nanome- 
ter scales, such as step bunching 1], meandering |2|, and 
pairing)^ on stepped surfaces under electromigration 
have been the focus of a large literature. In the usual 
model for step motion under electromigration^, mobile 
atoms between steps experience a constant drift force, 
while kinetic boundary conditions at steps are not af- 
fected by migration. It was first suggested in Ref.[5j that 
migration could lead to additional terms in step kinetic 
boundary conditions, and that these terms could modify 
the electromigration- induced drift of monolayer islands. 
Nevertheless, a microscopic explanation was lacking. 

From the work of Rous and Bly 6] it is known that, 
in the vicinity of steps, the backscattering of carriers 
strongly alters the migration force. In the present letter, 
we show that the boundary conditions used previously in 
the literature should be modified to account for the vari- 
ations of the migration force in the step region. To do 
so, we derive a step model from a phase field model with 
a non-constant migration force. The analysis of this step 
model reveals that the new contributions may produce all 
the instabilities known so far on vicinal surfaces. We fi- 
nally discuss the relevance of our results for experimental 
observations on various surfaces. 

We start with a microscopic description in terms of a 
phase field model Q. We use the simplest model, where 
the surface is described by two fields 4> and c, which re- 
spectively represent the normalized height of the surface, 
and the local concentration of mobile atoms. More so- 
phisticated models, which may account for arbitrary ki- 
netics at the steps are not considered here for the sake of 
simplicity Q. The dynamics of the phase field 4> is such 
that <p relaxes to a step and terrace structure. Terraces 
are wide regions where <fi is constant, and steps have a 
typical width W, as shown in Fig.l. The motion of the 
steps is driven by the departure from equilibrium in the 
step region, measured by the concentration variations: 
(c — c eq ), where c eq is the equilibrium concentration. The 



dynamics of is a relaxation with time-scale r^: 

r 5 t = -d^f + W 2 V 2 dp + \{c~c eq )d^g (1) 

The free energy density / is a periodic function of <p 
with minima for the values of <$> corresponding to the 
terraces. The coupling function g is also a function of if), 
with d$g > 0, and g^ — g+ — 1, where ± indicate the 
lower or upper side of the step respectively. The constant 
A controls the strength of the coupling between mobile 
atoms and steps. 

A global external force is present, related to the macro- 
scopic current in the bulk of the crystal. This current in- 
duces a local surface migration force M on mobile atoms. 
The force M depends on the surface local geometry via 
the phase field 4>. The resulting mass flux, oriented in 
the direction n e = M/|M|, reads: 
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where D is the local ^-dependent diffusion constant, and 
£ = fcsT/|M| is a 0-dependent lengthscale which char- 
acterizes the amplitude of the force. Local mass conser- 
vation then reads: 
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where h is the solid concentration, i.e. the number of 
solid atoms per unit area above a plane of arbitrary 
height parallel to the terraces. The variations of h are 
analogous to that of 0, and the jump of h across a step 
is constant h- — h+ = where fl is the atomic area. 
Writing Eq.®, we have discarded the tensorial coupling 
between the geometrical anisotropy of the step and the 
migration direction (i.e. the orientation of n e ), which 
could change the local migration direction in the vicinity 
of a step. We shall rather focus on another effect which 
is, in our opinion, the most important one: the varia- 
tion of the migration force amplitude in the step region. 
The phase field model presented here also neglects the 
dependence of the migration force on the adatom den- 
sity. Indeed, this effect should occur at high densities 
only 0, while the adatoms densities which are relevant 
for the experiments mentioned below are small. 
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FIG. 1: Phase field model with a non-constant migration 
force. The step region of width W is hatched. For the sake of 
simplicity, we have shown a ID cross-section, but the model 
in the text is 2D. Notations: c is the mobile atom concen- 
tration, (f> the phase field, D is the local diffusion constant, 
and |M| is the local amplitude of the migration force. Using 
suitable asymptotics, one obtains a step model as shown in 
the lower part of the figure, where the concentration may be 
discontinuous at the step (c+ 7^ c_), and <j> is now a Heaviside 
function which indicates the step position. 



We now derive a step model from the phase field model 
(11 IMI) . When step kinetics is fast, one may perform the 
thin interface asymptotics reported in Ref.0, Ei]- We 
shall assume that steps are slowly driven out of equi- 
librium, so that the smallest cut-off length l c related to 
the diffusion field on terraces is much larger than the 
step width W. We therefore define a small parameter 
e ~ W/£ c . We also take the coupling constant X ~ 1, 
and the departure from equilibrium (c — c eq ) ~ e. The 
latter condition indicates that step kinetics is fast enough 
to keep the concentration in the step region close to c eq . 
Furthermore, the relaxation dynamics of the phase field 
itself must be fast enough so that the dynamics does 
not affect the step and terrace structure, and we choose 
1> ~ e 2 . 

From these asymptotics, the dynamics of the concen- 
tration of atoms far from the steps reduce to Eq. © , with 
D — > Do, and without the term dth. 

On the macroscopic level, and in the limit where the 
adatom coverage is small, mass conservation at the step 
reads: 
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where £0 and Dq are respectively the constant value of 
£ and D on terraces far from steps. Performing the thin 
interface asymptotics, we obtain Eq.Q to leading order. 
But since c ~ c eq at the steps to leading order, the last 
term (c_ — c + )/£ s only appears to higher orders. 

The additional equations obtained from the thin inter- 
face asymptotics are the kinetic boundary conditions for 
the concentration at the steps 0: 



L> n.Vc± - 



D oc± ( - a V n , n c eq 

= v± I c± -c eq -f3— ±Q±yHj5) 



where c eq = c eq (l + Tn). These expressions involve some 
coefficients, the definition of which follows: 12JI 
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(7) 



where f m in is the value of / at its minima. Integrals are 
taken through one step in the normal direction z. The 
constants g± are the values of g on both sides of the steps. 

If we ignore the terms proportional to Q±, Eqs.JHI) are 
equivalent to the general linear kinetic boundary condi- 
tion of steps when step transparency and attachment- 
detachment on both sides of the steps are taken into 
account 7J. Using linearized non-equilibrium thermody- 
namics, Eqs.© could have been stated directly, with the 
new terms proportional to Q± [25J. Hence, Eqs.© are 
more general than the specific limit of the specific phase 
field model derived above. The relevance of our deriva- 
tion from a phase field model comes from the explicit 
link between the terms of Eqs.© and microscopic quan- 
tities. Indeed, from Eqs.0, Q+ or Q_ account for the 
variations of |M on the lower or upper side of the step 
respectively. Moreover, Q± > or Q± < respectively 
account for the a decrease or an increase of |M| in the 
step region. 

Instead of analyzing the consequences of the full 
boundary conditions 10, we shall now focus on a sim- 
plified limit. We define: 



(8) 



Q_ = / dx (I - V£) 



In order to grasp some of the main consequences of the 
new contributions, we shall focus on a simplified limit 
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where Q + = = Q/2. This may be obtained in the 
phase field model with an antisymmetric g and a sym- 
metric variation of l/£. This assumption amounts to 
considering that the main effect of the force variation in 
the step region is an average increase, or decrease of the 
amplitude of the force, respectively leading to a negative 
or positive Q. We also focus on a specific limit for step 
kinetics. We define: 



dx{D /D - 1) 



(9) 



Comparing the last term in the l.h.s. with the last term 
in the r.h.s. of Eqs.©, we infer that Q± will be relevant 
when it is of the order magnitude as Do jv± . We therefore 
take the limit where v± is large, which implies do small. 
We also take D to be symmetric across the step so that 
z/_i_ = V-. 

With these assumptions -symmetric l/£ and D and 
fast kinetics- the boundary conditions © read: 



c± = c eq (l + Tkt Q*/26) + PVn/tt 



where 



Q* = Q + d = dx 



Do 
D 



(10) 



(11) 



The expression indicates in a remarkably compact 
formulation, that a single quantity Q* accounts at the 
same time for the variations of the migration force, and 
for the variations of the diffusion constant in the step 
region. Eqs. QD} account for transparent steps when (3 is 
large, and non-transparent steps with fast attachment- 
detachment kinetics when /3 is small. 

We shall now point out how all known elementary in- 
stabilities of vicinal surfaces are obtained from the model 
with kinetic boundary conditions l|10|l . Let us analyze the 
stability of an array of initially straight and parallel steps 
separated by the same distance I. We assume that the 
direction of migration is orthogonal to the average step 
direction. We also consider the conserved regime where 
F = and 1/r — > 0. The mth step is perturbed by a 
small deviation 



Cm(x, t) = Cu4> exp[iu;t + irruj} + iqx] 



(12) 



where 0, and q are respectively the phase shift from step 
to step and the wavevector along the step. An instability 
is indicated by a positive 5Re[zu>]. Substituting Ea. (|12|) in 
the model equations, one finds the general expression of 
iuj as a function of 4> and q. 

In the case of in phase meandering (0 = 0), one finds 
in the long wavelength and weak electromigration limit 
((^) 2 ~7£o«l): 



iu ~ VLc eq D 



?V - T£q 4 



(13) 



Therefore, the train of steps is unstable when Q* /£ > 
0, and the most unstable wavelength is: A m = 
27r(2£ r£/Q*) 1 / 2 . 

In order to analyze step bunching, we shall introduce 
a repulsion between steps. If this interaction is of elastic 
or entropic origin, its free energy per unit step length is 
|ll| . The local equilibrium concentration is then 



A/t 



{1 + TK + 2A (£+ 3 -tZ 3 )) 



(14) 



where we have only considered the interaction between 
neighboring steps, and A — flA/hsT. For the sake of 
simplicity, we analyze the stability of the mode q = 0, 
which accounts for the bunching of straight steps. In 
the limit of weak electromigration, the linear stability 
analysis leads to: 

iu = (2Q*/£ l 2 )[l - cos(^)] - (12A/£ 5 )[1 - cos(</>)] 2 
nDc eq 1 + 2[1 - cos(4>)]nDf3/£ " ' 



Once again, an instability appears if > 0- The 

stability analysis thus indicates that a vicinal surface is 
simultaneously destabilized with respect to bunching and 
meandering for a downhill flux when migration or diffu- 
sion is weaker in the step region, and for an uphill mi- 
gration when migration or diffusion is enhanced in the 
vicinity of the steps. 

Let us now discuss some of the experimental results of 
the literature. On Si(lll), atoms drift along the current 
direction 01 • Three temperature regimes, denoted 1,11, 
and II, with increasing temperature, were found above 
the (7x7) — > (lxl) reconstruction transition tem- 
perature 0, 0| . In regimes I and III, bunching is ob- 
served for a downhill current during both growth and 
sublimation 1 13l| . In regime II, bunching is observed dur- 
ing growth for a downhill current, and during sublimation 
for an uphill current |l3| . From the step model without 
the Q contribution, it was concluded that regimes I and 
III correspond to opaque steps ((3, v + , and V- small), 
and regime II corresponds to fast step kinetics (f+, and 
i/_ large) 00. 

In the case of opaque steps, the Q terms are probably 
only a small contribution, unless a dramatic increase of 
migration is observed in a large region around the step. 
There is no evidence of such a strong effect. 

In the case of fast kinetics, the Q terms become impor- 
tant. From the study of step pairing Q, it was concluded 
that [3 is large, which means that steps are transparent. 
In Ref.Q, the occurrence of pairing was related to a small 
negative kinetic length do ~ — 0.13A. In the present anal- 
ysis, this result generalizes to Q* = Q + do « — 0.13A. 
We therefore obtain an alternative explanation for the 
pairing in the limit \Q\ ^> |do| of steps on Si(lll), based 
on Q < 0, i.e. on a stronger migration in the step regions. 
Quantitatively, we findQ « -0.13A. This may for exam- 
ple account for an increase of 1% of the migration force 
in a step region of width 13A. 
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Step meandering is also observed during sublimation 
with downhill migration in regime II ■ If Q < 0, as 
suggested above, then Q rather favors the stabilization 
of the meander for a downhill flux. Meandering therefore 
requires an extended study, which would also account for 
the sublimation rate, and which is beyond the scope of 
the present letter. A recent discussion can be found in 
Ref.0. 

Another system of interest is Si(100), which undergoes 
a dimer-row reconstruction, rotated of it/2 from one layer 
to the next one. We here assume that dimer rows are per- 
pendicular _L or parallel || to steps, and that the current 
is perpendicular to the average step orientation. It was 
shown in Ref. Q that the difference of diffusion constant 
(Dj_, _D||) perpendicularly to the step leads to pairing, 
and subsequent bunching of pairs for both directions 
of the electric current. The crucial quantity on terraces 
is the migration induced mass flux j ~ c eq D/^, which 
is proportional to the product of the diffusion constant 
with the migration force. Therefore, the phenomena that 
are attributed to a variation of D can also be the re- 
sult of a variation of £. Thus, pairing and bunching of 
pairs can also be the result variations of £. More pre- 
cisely, stronger migration along dimer rows reproduces 
pairing and bunching of pairs for both directions of the 
current, as observed in experiments [lfij . As a conclu- 
sion for Si(100), the variations of the force cannot be 
neglected a priori. Furthermore, following Ref.0], one 
could consider pairs as effective steps, with an internal 
region where diffusion and the migration force may vary. 
We could then follow the above mentioned analysis to 
determine the effective value of Q. This would help to 
analyze the pair bunching dynamics. 

We also expect the effect presented above to be rel- 
evant for the case of metals. Indeed, from microscopic 
models on metallic surfaces |f| , strong variations of the 
migration force in the vicinity of the steps were found. 
Rous and Bly 6] , indicate a decrease up to ~ 50% in a re- 
gion ~ 20A for Na. Using Eq. JHl with an approximate in- 
tegration of the results of Ref.jfj, we find Q « lOA. This 
contribution is two order of magnitude larger than that 
found above for Si(lll), and should therefore have dras- 
tic consequences on the dynamics of metal surfaces under 
electromigration. Experiments on Au surfaces |20( indi- 
cate that in presence of an electric current, the Ehrlich- 
Schwoebel (ES) effect (a decrease of interlayer mass 
transport related to a low value of v ) disappears, and 
a significant increase of interlayer mass transport is ob- 
served. The disappearance of the ES effect is probably 
caused by the increase of the temperature due to the joule 
heating effect. Indeed, it was shown in Ref . |2ll| . that the 
ES effect disappears at high temperature on Au. There- 
fore, kinetics is fast, and the tendency to bunching in- 
dicated in Ref.|20j could be caused by the non-constant 
migration scenario presented above. 

Step bunching was also observed on some other metal 



surfaces |22(. Moreover, recent experiments |23| have 
shown that GaN(OOOl) vicinal surfaces also exhibit step 
bunching for a downhill electric current during epitax- 
ial growth. The non-constant migration force scenario 
should also be considered in these cases. But the lack of 
quantitative understanding of the microscopic processes 
prevents a precise analysis. 

In conclusion, we have shown that variations of the mi- 
gration force in the vicinity of the steps may lead to step 
bunching, meandering, or pairing during electromigra- 
tion. This mechanism defines a novel scenario for sur- 
face instabilities, which may compete or combine with 
the other destabilizing mechanisms analyzed in the lit- 
erature. Finally, we have discussed the relevance of our 
results for some semiconductor and metal surfaces. 
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